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Introduction 

In  an  effort  to  treat  breast  cancer,  selective  estrogen  receptor  modulators  (SERMs) 
have  been  used  to  modulate  the  estrogen-signaling  pathway  with  mixed  results  [1].  A 
classic  example  of  a  SERM  is  tamoxifen.  When  used  as  a  therapeutic  for  a  limited  time, 
tamoxifen  is  effective  in  disrupting  the  estrogen-signaling  pathway.  Unfortunately,  with 
prolonged  use,  breast  tumor  cells  become  resistant  to  tamoxifen  and  are  able  to  use  the 
bioactivated  metabolite  of  tamoxifen  to  interact  with  co-activators  that  activate  the 
estrogen-signaling  pathway,  reversing  its  original  role  [1],  An  alternative  therapeutic 
approach  is  to  target  the  binding  site  of  the  co-activator  protein.  Recent  studies  have 
shown  that  some  small  molecules  may  bind  in  sites  (e.g.  co-activator  site)  other  than 
the  estradiol  binding  site  [2]  and  still  disrupt  the  estrogen-signaling  pathway.  By  binding 
in  the  co-activator  site  while  estradiol  is  bound  in  the  estrogen  receptor  (ER)  ligand 
binding  domain  (LBD),  these  small  molecules  act  as  co-activator  binding  inhibitors 
(CBIs)  because  the  co-activator  proteins  can  no  longer  bind;  thus,  gene  transcription  is 
inhibited.  Potentially,  these  CBIs  can  act  as  a  new  therapeutics  against  environmental 
or  natural  agonists  of  ERa.  Quantitative  structure-activity  relationship  (QSAR)  studies 
have  been  used  to  develop  therapeutics  that  will  compete  and  bind  in  the  estradiol 
binding  site  of  the  ERa  LBD  [3-5].  Because  these  studies  have  focused  on  the  estradiol 
binding  site,  new  potential  ER  disruptors  that  bind  in  the  co-activator  site  have  been 
missed.  Our  proposal  focuses  on  developing  a  new  computational  approach  to  predict 
therapeutically  useful  ERa  disruptors  by  investigating  CBIs  binding  to  the  co-activator 
site  in  conjunction  with  estrogenic  compounds  bound  in  the  estradiol  site 

Body 

Task  1 :  Computationally  predict  the  relative  binding  of  CBIs  in  the  co-activator  site  and 
the  CBI  binding  properties. 

Three  estrogenic  compounds  were  chosen  as  positive  controls,  estradiol  (E2),  1- 
methyl-2-amino-6-phenylimadazo[4,5b]pyridine  (PhIP),  and  diethylstilbestrol  (DES). 
These  three  compounds  have  been  shown  to  be  capable  of  activating  the  estrogen 
receptor  and  producing  an  estrogen  specific  response  (see  Appendix  A  and  more  [6]). 
Other  compounds  that  have  been  shown  to  inhibit  the  function  of  ERa  (Appendix  A) 
were  chosen  as  possible  CBI’s.  These  compounds  are  aromatic  heterocyclic  amines 
and  known  mutagens  [7,  8]. 

The  entire  ERa-LDB  was  used  in  the  docking  calculations  by  including  it  in 
all  the  atom  interaction  grids.  As  shown  in  appendix  A  and  figure  1,  E2  docking 
recovered  the  crystallographic  binding  mode  with  a  probability  of  binding  in  the  native 
binding  cavity  of  -99%.  DES  also  showed  a  very  high  probability  for  binding  in  the  E2 
site  (94%).  Striking,  however,  was  the  discovery  of  a  second  binding  site  in  the  interior 
of  the  protein  (Task  1.1).  This  site,  which  has  only  been  discussed  twice  before  in  the 
literature,  was  found  because  the  whole  protein  was  included  in  the  docking 
calculations.  The  paper  that  described  the  original  crystal  structure  defined  ths  region  as 
a  solvent  channel  [9].  A  second  mention  of  this  cavity  in  the  literature  was  an  attempt  to 
explain  unusual  kinetics  experiments  [10].  At  this  writing,  no  other  evidence  for  small 
molecule  binding  in  the  solvent  cavity  has  been  published.  E2  and  DES  showed  less 
than  1%  probability  of  binding  in  this  solvent  cavity. 
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In  our  ranking,  we  found  that  compounds,  which  inhibited  the  activation  of  the 
ER,  showed  a  preference  in  binding  to  the  alternative  site  (Figure  2  and  activiation 
figure)  (Task  1.2). 


Boltzman  Binding  Distributions 


Figure  1 :  The  Boltzman  binding  distribution  of  E2  and  other  heterocyclic  amines  on  the  ERa-LDB  at  the 
native  binding  cavity  and  the  back  cavity.  Details  of  the  calculations  and  compounds  are  described  in 
appendix  A. 

The  end  of  the  alternative  binding  site  adjoins  residues  important  for  binding  the 
co-activator  protein  (Figure  3).  We  hypothesized  that  compounds  bound  in  this 
alternate  cavity  could  affect  the  dynamics  of  residues  on  the  surface  of  the  protein. 
Molecular  dynamics  simulations  of  highest  ranking  CBI  (3M5-PhlP)  showed  that 
significant  distortion  of  the  protein  surface  was  possible  when  the  alternate  cavity 
contains  small  molecules  (Task  1.3).  In  our  simulations  of  3M5-PhlP  bound  to  ERa- 
LDB  we  noticed  that  Met  357  was  very  dynamic;  the  side  chain  would  move  significantly 
relative  to  surrounding  residues  during  the  2ns  simulation.  This  side  chain  motion 
affected  the  surface  of  the  protein,  especially  in  the  region  required  for  co-activator 
binding  (Figure  4).  Two  separate  simulations  had  similar  results.  Control  simulations 
of  E2  and  the  Tif-2  peptide  from  the  co-activator  protein  are  in  progress.  Comparison  of 
the  control  simulations  against  simulations  of  the  CBI  bound  in  the  alternate  cavity  will 
allow  for  better  quantification  of  the  surface  distortions. 
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Figure  2:  ERa-LDB  shown  with  the  backbone  in  ribbons  and  binding  cavities  in  surface  representation 
colored  by  residue  type  (red,  negative  charge;  blue,  positive  charge;  green,  polar;  white  non-polar).  The 
protein  is  rotated  180°  to  show  the  opposite  positions  of  the  cavities  on  the  surface.  E2  is  shown  in  stick 
representation  inside  the  crystallographic  binding  site,  while  the  PhIP  congener  3-methyl-5phenyl-PhlP  is 
shown  in  its  predicted  “alternate  binding  cavity”.  The  blue  surface  on  the  protein  at  the  left  defines  the 
region  of  residues  critical  for  binding  the  co-activator  protein. 
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Figure  3:  The  ERa-LDB  with  a  bound  co-activator  peptide  NR2  from  TIF2.  The  ER  protein  backbone  is 
colored  in  cyan  and  shown  in  cartoon  with  cylinders  and  arrows  signifying  a-helices  and  (3-strands, 
respectively.  The  peptide  is  drawn  in  stick  representation.  Surface  residues  of  the  ERa  protein  that 
interact  with  the  peptide  are  shown  as  surfaces  colored  by  residue  type  as  described  in  figure  2.  Arrows 
show  positions  of  the  charge  clamp  residues,  Asp  542  and  Lys  362. 
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Figure  4:  The  ERa-LDB  before  molecular  dynamics  and  after  2ns  of  simulation  time.  The  protein 
backbone  is  shown  in  ribbon  representation  and  colored  in  cyan.  Methionine  357  is  shown  in  stick 
representation.  Surface  residues  that  interact  with  the  helical  peptide  of  the  co-activator  protein  are  in 
blue.  Residues  that  become  accessible  to  solvent  after  simulation  are  shown  in  magenta.  The  difference 
in  exposed  surface  area  of  the  protein  is  about  200  A2.  Before  the  simulation  is  started  the  sidechain  of 
Methionine  357  points  in  towards  the  alternate  cavity  were  3M5-PhlP  is  present.  Early  in  the  simulation, 
the  sidechain  switched  positions  to  point  out  towards  the  solvent  accessible  surface.  This  disrupts  at 
least  one  small  hydrophobic  pocket  that  is  required  for  optimal  co-activator  binding. 

The  resulting  combination  of  CBIs  and  estrogenic  compounds  are  currently  being 
re-ranked  based  on  MD  simulations.  We  have  also  received  a  40,000  compound  library 
to  screen  using  our  methods  developed  in  this  project  and  this  work  is  continuing.  (Task 

1.4) 

Task  2:  Biologically  assay  the  effect  of  the  CBI  on  the  estrogen  receptor  activation. 

The  ability  of  potentially  estrogenic  compounds  to  activate  the  estrogen  receptor 
was  assessed  using  a  standard  luciferase  assay.  This  assay,  a  discussion  of  the  cells 
used  for  this  experiment,  and  results  of  the  effects  of  several  heterocyclic  amines  on  the 
ER  activation  are  described  in  Appendix  A  (Tasks  2.1  and  2.2). 

In  addition,  using  this  assay,  we  were  able  to  assess  the  effects  of  small 
molecule  structure  on  ER  activation.  Using  a  series  of  heterocyclic  amine  compounds 
that  differed  only  slightly  in  structure,  we  were  able  to  show  that  PhIP  was  optimally 
configured  to  activate  the  ER  (figure  5). 


Control  Estradiol  Genistein  PhIP  N-OH-PhIP  3-Methyl  3-Methyl-5-  IQ  MelQx  IFP  TMIP 

PhIP  Phenyl 
PhIP 

Figure  5:  Effect  of  heterocyclic  amines  on  ER  activation.  Concentrations  of  all  HAs  tested  was  4x10'7M. 

E2  concentrations  was  10"9  M.  A  discussion  of  the  relevance  of  including  genistein  in  this  experiment  is 
provided  in  Appendix  A. 

Competition  assays  were  used  to  determine  the  effect  of  simultaneous  incubation 
of  HAs  (4x10"7M)  and  10'9  M  E2  on  ERa  activation  (Fig.  6).  Treating  MCF-7  cells  with 
PhIP  and  E2  produces  a  response  that  is  7%  less  than  E2  alone.  Ar-Hydroxy-PhIP 
inhibits  the  E2-mediated  activation  of  ERa  by  30%,  and  the  other  HAs,  with  the 
exception  of  TMIP,  reduce  E2-initiated  Er  activation  to  the  levels  of  the  control  cells. 
Based  on  the  docking  data,  we  find  that  the  compounds  that  are  most  strongly  predicted 
to  bind  to  the  alternative  binding  cavity  on  the  ERa  protein,  are  the  best  inhibitors  (Task 
2.3) 
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Figure  6.  Luciferase  Activation  Based  Competition  assays,  HAs  and  E2.  Concentrations  of  all  HAs  tested 
were  4x1 0‘7M.  E2  concentrations  were  1 CT9  M.  Total  activation  relative  to  E2  on  the  y-axis  vs  HAs  on  the 
x-axis. 


We  wanted  to  know  how  some  HAs  were  inhibiting  ERa  activation  while  PhIP  did 
not.  A  simple  competition  assay  was  used  in  which  E2  plus  an  HA  were  added 
simultaneously  to  the  incubated  cells.  The  resulting  data  are  shown  in  Figure  6. 
Estrogen  overcame  the  inhibitory  effect  of  2NOH-PhlP  and  TMP  while  in  enhancing  the 
activating  properties  of  PhIP.  Other  HAs  were  close  to  the  control  activity.  These  data 
lead  us  to  believe  that  an  alternative  binding  site  is  present  at  the  level  of  the  ERa-LDB. 
It  is  possible  that  other  steps  in  the  ER  signaling  pathway  are  the  actual  targets  of  the 
inhibitory  HAs,  NMR  experiments  are  planned  to  confirm  or  deny  binding  to  the  ERa 
protein. 
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Task  3:  Spectroscopically  measure  the  binding  of  the  CBI  to  the  ER  LBD 

Full  length,  recombinant  human  estrogen  receptor  a  (hERa)  is  commercially 
available  through  Sigma  Aldrich,  Inc.  It  is  supplied  as  1 .5  pM  purified  protein  (750 
pmoles  in  500  pi  of  50  mM  tris-HCI,  pH  8.0,  500  mM  KCI,  2  mM  EDTA,  1  mM  sodium 
orthovanadate  and  10%  glycerol).  Recombinant  hERa  is  produced  as  an  active,  soluble 
66.4  kDa  protein  by  a  baculovirus  expression  system,  and  thus  possesses  post- 
translational  modification.  For  NMR  binding  studies,  it  was  first  necessary  to  determine 
alternative  buffer  conditions  in  which  the  ER  protein  remained  stable.  Both  tris  and 
glycerol  contribute  large  proton  signals  in  the  NMR  data  that  interfere  with  the  detection 
of  submicromolar  quantities  of  sample,  in  addition,  high  salt  concentrations  can  result  in 
protein  precipitation  when  small  volumes  of  compounds,  such  as  PhIP  and  E2,  which 
are  only  soluble  in  organic  solvents  such  as  DMSO,  are  added.  After  numerous  trials, 
we  had  determined  that  the  protein  was  able  to  remain  stable  for  a  few  days  after  being 
exchanged  after  extensive  dialysis  (Spectrapor  7,  MWCO  1000  dialysis  tubing)  into 
10mM  Na2HP04,  pH  8.0,  120mM  NaCI,  2.7mM  KCI  (PBS)  buffer  (Task  3.1). 

Various  NMR  spectroscopic  techniques  can  be  used  to  screen  for  interactions 
between  proteins  and  ligands  (reviewed  in  [11,  12].  Transferred  NOESY  (trNOESY) 
experiments  are  routinely  used  to  detect  ligand  binding  to  a  target  protein  under 
conditions  of  fast  exchange  (ligands  that  bind  with  mM  to  mM  dissociation  constants) 
[13].  The  intensity  of  each  intra-ligand  NOE  crosspeak  is  governed  by  the  population- 
weighted  cross-relaxation  rate.  A  strong  negative  NOE  crosspeak  is  observed  for 
binders  (black  peaks),  as  opposed  to  weakly  positive  (red  peaks)  or  zero  NOE 
crosspeaks  for  non-binders  or  in  the  absence  of  protein.  Thus,  the  sign  flip  of  the  NOE 
cross  peak  between  the  free  versus  bound  states  acts  as  a  simple  binary  filter  to 
distinguish  binders  from  non-binders. 
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mM  estradiol  dissolved  in  deuterated  DMSO.(c)  The  expanded  region  of  a  900  ms  mixing  time  NOESY 
spectrum  of  5.95  mM  PhIP  in  DMSO  exhibits  weak  positive  NOE  crosspeaks  (boxed  in  red),  (d)  The 
expanded  region  of  a  300  ms  mixing  time  NOESY  spectrum  of  3  pM  ER  and  120  pM  PhIP  indicates  that 
PhIP  binds  as  evidenced  by  the  strong,  negative  crosspeaks  (boxed  in  black). Concentration  of  PhIP  is 
50X  greater  in  (c)  than  in  (d).  Spectrum  in  (d)  is  plotted  at  2  X  lower  intensity  level  relative  to  spectrum  in 
(c)  for  presentation  purposes,  (e)  The  expanded  300  ms  mixing  time  NOESY  spectrum  after  addition  of 
1 .5  pM  estradiol  to  the  ERa+PhIP  sample  shows  that  PhIP  does  not  bind  (red  peaks).  Spectrum  (e)  is 
plotted  at  2X  lower  intensity  level  relative  to  spectrum  in  (d).  This  result  suggests  that  PhIP  has  been 
displaced  from  binding  to  ERa  by  estradiol,  and  that  both  ligands  may  be  binding  in  the  same  site.  Note 
that  the  chemical  shifts  of  PhIP  change  depending  on  solvent  (DMSO  in  spectrum  (c)  versus  aqueous 
buffer  in  spectra  (d)  and  (e)). 


We  observed  a  sign  flip  in  NOE  crosspeaks  of  PhIP  from  positive  (red)  to 
negative  (black)  in  the  trNOESY  spectra  when  PhIP  was  added  to  ERa,  confirming 
experimentally  that  PhIP  does  indeed  bind  to  ER  (Figure  7c  and  d)  (Task  3.2).  Since 
E2  is  known  to  bind  with  high  affinity  to  ERa  (Kd  ~  10'8-1CT9  M),  we  could  not  use  the 
trNOESY  method  to  detect  its  binding  to  ERa.  However,  since  the  E2’s  binding  site  on 
ERa  is  known,  it  is  possible  to  use  this  information  to  design  a  trNOESY  competition 
experiment  that  can  provide  information  about  where  PhIP  is  binding.  If  PhIP  is  binding 
in  a  different  site  than  estradiol,  than  we  should  be  able  to  still  observe  negative  (black) 
peaks  because  PhIP  and  estradiol  are  binding  in  two  different  sites,  and  estradiol  is  not 
displacing  PhIP  binding.  In  contrast,  if  PhIP  is  binding  in  the  same  site  as  estradiol, 
then  adding  estradiol  to  the  PhIP/ERa  mixture  should  competitively  displace  it  since 
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estradiol  binds  much  more  strongly  to  ERa.  In  this  case,  we  would  observe  another 
sign  flip  of  the  PhIP  peaks  from  black  (binding)  to  red  (not  binding).  Addition  of  estradiol 
to  the  PhIP/ERa  mixture  resulted  in  PhIP  cross  peaks  flipping  sign  from  black  to  red 
(Figure  7d  and  e)  and  thus  suggesting  that  both  PhIP  and  estradiol  are  binding  in  the 
same  site  on  Era  (Task  3.3).  This  result  is  consistent  with  the  results  shown  by  us  and 
others  that  PhIP  exposure  increases  MCF-7  cell  proliferation  and  ERa  activation.  It  also 
agrees  with  our  computational  model  and  with  the  idea  that  PhIP  directly  binding  to  ERa 
is  responsible  for  the  effects  observed.  However,  we  can  not  completely  rule  out  the 
possibility  that  PhIP  is  binding  in  another  site,  and  that  estradiol  prevents  PhIP  from 
binding  because  of  a  conformation  change  in  the  protein. 

Thus  in  order  to  further  clarify  the  identity  of  the  ligand  binding  site,  we  initiated 
the  expression  of  15N-labeled  ERa  ligand  binding  domain  (LBD)  (Task  3.4).  Dr.  Myles 
Brown  (Harvard  University)  has  generously  provided  us  with  the  expression  vector  for 
the  hERa-LBD).  The  GST-LBD  fusion  protein  exhibits  binding  affinity  to  estradiol  (KD  ~ 
0.1  to  3.3  nM)  that  is  comparable  with  reported  values  measured  for  wild-type  MCF-7 
ERa  expressed  in  vivo  [14]. 

The  hERa-LBD  is  expressed  as  a  GST  fusion  protein  in  BL21  (DE3)  cells.  Cell 
growth  is  carried  out  for  16  hours  at  37  SC  with  shaking  at  300  rpm  (INOVA  shaker) 
using  either  an  autoinduction  method,  as  introduced  by  Studier  [15]  or  with  a  modified 
version  of  the  autoinduction  method  described  specifically  for  NMR  studies  [15,  16]. 
Cells  are  harvested  by  centrifugation  and  lysed  by  sonication  on  ice  in  a  buffer 
containing  50mM  Tris-CI,  pH  7.5,  lOOmM  NaCI,  20mM  [5-  mercaptoethanol,  0.5%  NP- 
40,  followed  by  centrifigation  at  30,000g  at  45C  for  30min  to  remove  cell  debris. 

Although  the  expression  yields  are  very  high  for  the  GST-ERa-LBD  protein,  as 
can  be  seen  in  Figure  8,  the  protein  is  insoluble.  Currently,  we  are  in  the  process  of 
modifying  the  expression  protocol  (addition  of  glycerol,  tween  and  other  detergents,  as 
well  as  lower  temperature)  to  increase  the  solubility  of  the  protein.  Once  this  is 
accomplished,  the  resulting  supernatant  will  then  be  diluted  and  bound  to  a  glutathione 
sepharose  4B  column  (Amersham)  [17].  After  several  washes  with  PBS  buffer,  the 
hERa-LBD  protein  will  be  cleaved  from  GST  and  eluted  from  the  sepharose-GST 
column  by  treatment  with  thrombin  (Amersham).  Yields  will  be  determined  by  the 
Bradford  protein  assay  (BioRad).  The  purified  protein  will  then  be  concentrated  into 
PBS  buffer,  pH  7.4  by  centrifugation  using  either  a  Centricon  YM-3  or  Ultra-4 
concentrators  (Amicon)  [16]. 
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Figure  8:  Polyacrylamide  gel  electrophoresis  (4-12%)  results  showing  high  expression  levels  of  GST 
ERa-LBD.  (a)  Most  of  the  protein  is  in  the  insoluble  fraction  (lane  1 :  GST-EFta-LBD  uninduced;  lane  2: 
induced  1  hr.  @  37C;  lane  3:  induced  2hrs.  @  37C;  lane  4.  clear  lysate  (soluble  fraction);  lane  5.  insoluble 
fraction,  (b)  ERa  expression  in  autoinducing  minimal  media,  (lane  1 :  molecular  weight  markers;  lane  2: 
Grew  ~16hrs,  37C,  300rpm,  cell  lysate;  lane  3: 10  times  dilution). 

The  availability  of  15N-labeled  ERa-LBD  protein  will  allow  us  to  carry  out 
Structure  Activity  Relationship  by  NMR  (SAR-by-NMR)  studies  [18]  that  will  clearly 
define  the  binding  site  and  orientation  of  not  only  PhIP,  but  other  CBIs  and  HAs. 

Task  4:  Submit  manuscript  for  publication. 

A  large  portion  of  the  work  detailed  above  in  tasks  1-3  has  been  submitted  and 
accepted  for  publication  in  the  American  Chemical  Society  journal,  Chemical  Research 
in  Toxicology.  At  the  time  of  this  report,  the  projected  publication  date  is  October  5th 
2005.  We  anticipate  submission  of  a  second  manuscript  describing  the  characterization 
of  the  co-activator  binding  site  interactions  with  ligands  bound  in  the  alternative  cavity 
through  15N  NMR  and  molecular  dynamics  simulations. 
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Key  Research  Accomplishments 

Posters  and  Accepted  Abstracts: 

•  19th  National  Meeting  of  the  Protein  Society  San  Diego  CA,  August  2004. 

•  K.S.  Kulp,  B.J.  Bennion,  J.L.  Montgomery,  F.L.  Lightstone,  M.G.  Knize,  J.S. 
Felton  and  L.M.  Bennett  “Dietary  constituents  affect  receptor  activation  and  cell 
proliferation  in  human  breast  cancer  cells”.  American  Association  of  Cancer 
Research  Annual  Meeting,  March  27-31, 2004,  Orlando  Florida. 

•  B.J.  Bennion,  K.S.  Kulp,  J.L.  Montgomery,  F.L.  Lightstone,  M.G.  Knize,  J.S. 
Felton  and  L.M.  Bennett  “Dietary  constituents  affect  receptor  activation  and  cell 
proliferation  in  human  breast  cancer  cells”.  American  Chemical  Society  Annual 
Meeting,  March  27-31 , 2004,  Anaheim  California. 

•  DOD  BCRP  Meeting  “Era  of  Hope”  Philadelphia  PA,  June  2005. 

•  36th  Meeting  of  The  Environmental  Mutagen  Society  San  Francisco  CA, 
September  2005. 

Manuscripts: 

•  “PhIP  Carcinogenicity  in  Breast  Cancer:Computational  and  Experimental 
Evidence  for  Competitive  Interactions  with  Human  Estrogen  Receptor  “  in  press 
Chemical  Research  in  Toxicology,  Oct  2005. 

Reportable  Outcomes 

Presentations: 

•  Biosciences  Directorate  Symposium  LLNL,  November  2004. 

•  Biosciences  Directorate  Postdoctoral  Symposium  LLNL,  July  2005. 

Invited  Talks-Lectures: 

•  University  of  California-Davis  Cancer  Center,  Sacremento  CA,  January  2005. 

•  Chemistry  Department  University  of  the  Pacific,  Stockton  CA,  February  2005. 

•  Bio-engineering  Department  University  of  the  Pacific,  Stockton  CA,  March  2005. 

•  Edward  Teller  Education  Center,  Lawrence  Livermore  Natl.  Lab.,  July  2005. 

NIH  R01  grant  entitled  “Dietary  exposure  to  multiple  heterocyclic  amines  may  cause 
fewer  breast  tumors  than  exposure  to  single  carcinogens”  to  be  submitted  October  1 , 
2005. 

Conclusions 

Computer  aided  docking  discovered  a  food  mutagen  that  could  bind  in  the  native 
binding  cavity  of  ERa-LDB  and  possibly  activate  the  estrogen-signaling  pathway.  The 
computational  work  also  predicted  a  possible  secondary  binding  site  for  small 
molecules.  In  deed,  we  showed  that  disruption  of  the  surface  of  the  protein  was 
possible  when  this  alternate  binding  site  is  occupied  by  small  heterocyclic  amines. 
Future  work  in  this  area  would  include  a  computationally  based  screening  of  a  larger 
database  of  commercially  available  compounds  that  are  similar  size  and  structure  to  the 
heterocyclic  amines  that  antagonized  ERa  gene  expression. 
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PhIP  was  also  shown  to  increase  proliferation  of  MCF-7  breast  cancer  cells  as 
well  as  activate  ERa  dependent  gene  expression.  Other  PhIP  related  heterocyclic 
amines  decreased  breast  cancer  cell  growth  and  blocked  ERa  activation  by  estradiol. 
These  findings  corroborate  the  predictions  made  in  task  1 ,  i.e.  PhIP  was  a  weak  agonist 
and  that  the  other  heterocyclic  amines  were  strong  antagonists  to  ERa  dependent  gene 
expression. 

NMR  experiments  showed  that  PhIP  could  bind  to  human  ERa.  E2  was  also 
shown  to  displace  bound  PhIP  from  human  Era.  These  two  observations  substantiated 
computational  predictions  in  task  1 ,  i.e.  P  hIP  was  predicted  to  bind  to  the  native  binding 
cavity  at  a  reduced  affinity  compared  to  E2.  In  addition,  we  also  developed  expression 
techniques  for  15N-labeled  human  ERa-LDB  protein  to  accurately  define  the  binding 
location  of  small  molecules.  This  research  has  resulted  in  a  submitted  manuscript  that 
was  accepted  for  publication  in  Chemical  Research  in  Toxicology  in  2005. 

In  conclusion,  we  have  shown  that  computational  docking  has  predictive  value,  which 
can  open  new  lines  of  research;  moreover  a  possible  new  link  between  food  mutagens 
and  breast  cancer  was  discovered. 
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Abstract 


Many  carcinogens  have  been  shown  to  cause  tissue  specific  tumors  in  animal  models. 
The  mechanism  for  this  specificity  has  not  been  fully  elucidated  and  is  usually  attributed 
to  differences  in  organ  metabolism.  For  heterocyclic  amines,  potent  carcinogens  that 
are  formed  in  well-done  meat,  the  ability  to  either  bind  to  the  estrogen  receptor  and 
activate  or  inhibit  an  estrogenic  response  will  have  a  major  impact  on  carcinogenicity. 
Here  we  describe  our  work  with  the  human  estrogen  receptor  alpha  (ERa)  and  the 
mutagenic/carcinogenic  heterocyclic  amines  PhIP,  MelQx,  IFP,  and  the  hydroxylated 
metabolite  of  PhIP,  A^-hydroxy-PhlP.  We  demonstrate  that  PhIP  binds  with  the  ligand 
binding  domain  (LBD)  both  by  computational  docking  and  NMR  analysis.  This  binding 
competes  with  estradiol  (E2)  in  the  native  E2  binding  cavity  of  the  receptor.  In  in  vitro 
assays,  we  find  that  PhIP,  in  contrast  to  the  other  heterocyclic  amines,  increases  cell- 
proliferation  in  MCF-7  human  breast  cancer  cells  and  activates  the  ERa  receptor.  We 
also  find  that  other  heterocyclic  amines  and  A^-hydroxy-PhlP  inhibit  ERa  activation.  We 
propose  that  the  mechanism  for  the  tissue  specific  carcinogenicity  seen  in  the  rat  breast 
tumors  and  the  presumptive  human  breast  cancer  associated  with  the  consumption  of 
well-done  meat  maybe  mediated  by  this  receptor  activation. 


Introduction 

Changes  in  breast  cancer  incidence  among  immigrant  populations  suggest  that 
lifestyle  factors,  including  diet,  may  be  an  important  cause  of  the  disease  and  a  potent 
clue  for  treatment.  One  dietary  modification  that  is  frequently  found  among  immigrants 
from  eastern  to  western  countries  is  an  increase  in  the  consumption  of  cooked  muscle 
meats.  Well-done  cooked  muscle  meats  are  known  to  contain  potent  mutagens  and 
mammary  carcinogens  (e.g.  in  rodents)  belonging  to  the  heterocyclic  amine  (HA)  class 
of  chemical  compounds  (Felton  et  al.,  2004).  Three  of  the  heterocyclic  amines 
commonly  found  in  meats  cooked  under  household  conditions  are:  2-amino-1-methyl-6- 
phenylimidazo[4,5-jb]pyridine  (PhIP),  2-amino-3,8-dimethylimidazo[4,5-/]quinoxaline 
(MelQx),  and  2-amino-1 ,6-dimethylfuro[3,2-e]imidazo[4,5-b]pyridine,  (IFP).  The  relative 
amounts  of  these  compounds  formed  during  cooking  depend  on  both  meat  type  and 
cooking  conditions  (Knize  et  al.,  1998).  PhIP  is  frequently  the  most  mass  abundant 
heterocyclic  amine  produced  during  the  cooking  of  beef,  pork,  and  chicken,  (Keating  et 
al.,  2000;  Knize  et  al.,  1998;  Norrish  et  al.,  1999;  Pais  et  al.,  2000;  Sinha  et  al.,  1995; 
Skog  et  al.,  1997;  Wakabayashi  et  al.,  1992).  Humans  are  routinely  exposed  to  varying 
amounts  of  these  food-derived  compounds,  and  there  are  studies  supporting  their  role 
in  human  carcinogenesis  (Knize  and  Felton,  2005). 

When  given  in  the  diets  of  rats,  PhIP  has  been  shown  to  cause  the  formation  of 
hormone-dependent  mammary  tumors  (Ito  et  al.,  1991).  A  powerful  liver  carcinogen, 
MelQx  also  causes  breast  tumors  in  Sprague-Dawley  rats  but  is  a  much  less  potent 
mammary  carcinogen  than  PhIP  (Snyderwine,  2002;  Wakabayashi  et  al.,  1992).  IFP 
has  been  shown  to  be  a  potent  mutagen,  but  has  not  yet  been  tested  for  carcinogenicity 
in  an  animal  model.  Although  the  mechanism  of  carcinogenesis  for  these  compounds 


has  not  been  fully  elucidated,  metabolic  activation  and  subsequent  formation  of  DNA 
adducts  is  believed  to  be  critical  (Snyderwine,  2002;  Snyderwine  et  al.,  2003).  The 
metabolic  activation  of  PhIP,  a  two-phase  process,  is  representative  of  the  activation 
pathways  of  the  other  HAs.  During  Phase  I  metabolism,  PhIP  is  oxidized  via  cytochrome 
P4501A2  (CYP1 A2)  to  a  hydroxylated  intermediate,  2-hydroxyamino-1-methyl-6- 
phenylimidazo  [4,5 -b\  pyridine  (A/^-hydroxy-PhlP).  A/^-Hydroxy-PhlP  is  then  converted  to 
a  more  biologically  reactive  form  via  Phase  II  metabolizing  enzymes,  primarily  the 
acetyltransferases  or  sulfotransferases.  This  esterification  generates  electrophilic  O- 
sulfonyl  and  O-acetyl  esters,  which  have  the  capacity  to  bind  DNA  and  cellular  proteins. 

In  addition  to  diet,  hormonal  factors  also  play  a  role  in  mammary  gland 
tumorigenesis.  Estrogen  receptor  a  (ERa),  a  member  of  the  nuclear  receptor  family, 
regulates  estrogen-responsive  gene  expression  upon  ligand  binding.  Structurally,  the 
ERa  protein  contains  four  major  domains:  activation,  DNA  binding,  hinge,  and  ligand 
binding  domains  (LBD).  Crystal  structures  of  LBD-ERa  with  agonist  and  antagonists 
bound  have  been  published  (Eiler  et  al.,  2001 ;  Henke  et  al.,  2002;  Pike  et  al.,  1 999; 

Pike  et  al.,  2001 ;  Renaud  et  al.,  2003;  Ruff,  1999;  Shiau  et  al.,  1998;  Shiau  et  al.,  2002; 
Warnmark  et  al.,  2001).  Ligand  dependent  ERa  activation  begins  with  17-(3-estradiol 
(E2)  binding  the  ERa  monomer  and  releasing  protein  chaperones  from  the  receptor 
(Pratt  and  Toft,  1997).  With  E2  bound,  ERa  monomers  then  dimerize  (Kumar  and 
Chambon,  1 988)  and  interact  with  co-activator  and  co-repressor  proteins  (Moras  and 
Gronemeyer,  1998).  Transcription  begins  when  all  factors  have  been  recruited  and  ERa 
binds  to  the  estrogen  response  element  (ERE)  of  the  specific  gene. 


Because  exposure  to  PhIP  alone  is  sufficient  to  cause  tumors  in  rodent  models, 
PhIP  may  be  important  in  mammary  tumor  development  beyond  general  DNA  adduction 
and  tumor  initiation.  Exposure  to  PhIP  has  been  shown  to  increase  cell  proliferation  in 
mammary  gland  terminal  end  buds,  potential  sites  of  tumor  development,  suggesting 
that  PhIP  promotes  tumorigenesis  by  causing  the  further  replication  of  PhIP-DNA 
adducts  in  target  cells  (Snyderwine,  1999).  This  potential  role  in  promotion  is 
strengthened  by  the  findings  of  Pfau  et  al.  and  Gooderham  et  al.  demonstrating  that 
PhIP  exposure  increases  human  breast  cancer  cell  proliferation  and  activates  ERa 
(Gooderham  et  al.,  2002;  Pfau  et  al.,  2000)  at  concentrations  that  are  less  than  those 
required  for  mutagenic  activity.  Gooderham  and  coworkers  have  also  shown  that  the 
pure  anti-estrogen  ICI  182,780  abolishes  PhIP  transcriptional  activity,  further  suggesting 
that  PhIP  binds  specifically  to  ERa  (Lauber  et  al.,  2004).  However,  in  studies  on 
ovarectomized  Sprague-Dawley  rats,  PhIP  did  not  produce  any  estrogenic  response 
(Kawamori  et  al.,  2001),  suggesting  that  PhIP  is  not  a  strong  estrogen  mimic  in  this 
assay  system.  Other  studies  propose  that  PhIP  promotes  carcinogenesis  by  retarding 
mammary  gland  development  (Snyderwine  et  al.,  1998),  increasing  serum  prolactin 
levels  (Venugopal  et  al.,  1999b)  and  inhibiting  apoptosis  (Venugopal  et  al.,  1999a). 

Here,  we  are  the  first  to  report  that  (1)  PhIP  binds  directly  to  the  ERa-ligand 
binding  domain  and  competes  with  estradiol  in  the  native  binding  cavity,  and  (2)  the 
PhIP  phase  I  metabolite,  A/^-hydroxy-PhlP,  and  two  heterocyclic  amines,  MelQx  and 
IFP,  inhibit  ERa  activation.  We  also  confirm  the  proliferative  effect  of  PhIP  on  MCF-7 
cells  and  the  estrogen  receptor  activation  shown  by  Lauber  et  al.  (Lauber  et  al.,  2004). 


Taken  together,  these  results  suggest  that  PhIP  activates  the  estrogenic  response  via 
direct  binding  to  ERa,  thereby  explaining  its  tumor  specificity. 

Methods 

Ab  initio  calculations 

Each  ligand  (E2,  PhIP,  A^-hydroxy-PhlP,  IFP,  and  MelQx)  was  geometry 
optimized  at  the  RHF  6-31 G*  level  with  Gaussian  98,  revision  A.1 1 .4  (Frisch  M.  J.  et 
al.,  2002)  (Fig.  1  A).  Default  settings  for  convergence  criteria  were  satisfied  for  each 
ligand.  Atomic  partial  charges  for  the  ligands  were  determined  by  electrostatic  potential 
fitting  (Wang  et  al.,  2000). 

Homology  Modeling 

A  high-resolution  crystal  structure  of  the  ERa  ligand  binding  domain  (1GWR, 
(Warnmark  et  al.,  2002))  was  chosen  from  the  Protein  Data  Bank  to  model  protein- 
ligand  interactions  in  docking.  There  were  residues  in  two  loops  and  other  atoms  in  the 
structure  where  electron  density  was  not  present  in  chain  A  (Warnmark  et  al.,  2002). 
These  loops  and  missing  atoms  were  modeled  by  aligning  1GWR  to  another  ER 
structure  (PDB  ID  1G50),  using  an  in-house  alignment  program,  LGA  (Zemla,  2003). 
The  primary  sequences  of  the  modeled  loops  were  identical  between  structures. 
Appropriate  hydrogen  atoms  were  added  to  the  structure  using  the  CFIARMM  27  force 
field  (MacKerell  et  al.,  1998). 


Docking 


The  estrogen  receptor  protein  (1GWR  homology  model)  and  ligands  (E2,  Phi P, 
/'/-hydroxy-PhIP,  and  MelQx)  were  prepared  for  the  docking  studies  as  described  in  the 
Autodock  manual  (Morris  et  al.,  1998).  The  original  source  code  was  modified  to  allow  a 
grid  cube  size  of  128  points  with  a  0.375  A  grid  spacing.  The  resulting  grid  covered 
-80%  of  the  ligand  binding  domain.  Using  the  lamarkian  genetic  algorithm  option,  the 
total  number  of  docking  steps  was  set  to  200  and  cluster  analysis  was  activated.  Each 
ligand  pose  was  then  grouped  according  to  its  position  on  the  protein.  A  Boltzmann 
probability,  Ki,  was  then  derived  from  the  calculated  free  energy  AG|,  gas  constant  R, 
and  temperature,  T (30°  C)  of  each  pose,  /'. 


K.  =  <%AG‘  'RT) 


i  -  Equ.  1 

Docking  poses  in  each  grouping,  s,  were  then  weighted  according  to  Equ.  2. 


Equ.  2 


Here,  p,  is  the  weighted  probability  of  finding  pose,  /,  in  grouping,  s.  These  probabilities 
were  then  summed  for  each  category. 


Transferred  Nuclear  Overhauser  Effect  Spectroscopy  (trNOESY) 

All  spectra  were  recorded  at  30°  C  using  a  Varian  Inova  600  MHz  spectrometer.  ^-1 D 
and  phase  sensitive  2D-nuclear  Overhauser  effect  spectroscopy  (NOESY)  experiments 
(900  ms  mixing  time)  were  carried  out  on  1 .0  mg  of  PhIP  (5.94  mM)  or  27  mg  of  E2 


(132.1  mM)  dissolved  in  750  pL  deuterated  DMSO)  in  the  absence  of  ERa  protein  (Fig. 

1 B,  C).  Human,  recombinant  ERa  protein  was  purchased  from  Sigma-Aldrich  as  750 
pmols  of  purified  protein  in  500  pL  50  mM  Tris-HCL,  pH  8.0,  500  mM  KCL,  2  mM  DTT, 
ImM  EDTA,  1  mM  sodium  orthovanadate  and  10%  gylcerol.  We  first  determined 
alternative  buffer  conditions  in  which  the  ERa  protein  remained  stable.  Both  Tris  and 
glycerol  contribute  large  proton  signals  in  the  NMR  data  that  interfere  with  the  detection 
of  micromolar  quantities  of  PhIP  and  E2.  In  addition,  high  salt  concentrations  are  not 
compatible  with  adding  small  volumes  of  PhIP  and  E2,  which  are  only  soluble  in  organic 
solvents  such  as  DMSO,  and  results  in  precipitation  of  the  protein.  We  determined  that 
ERa  remained  stable  for  several  days  after  being  exchanged  (using  dialysis, 

MWCO=1 0,000)  into  50  mM  Na2P04,  pH  8.0,  200  mM  KCL,  1%  p-mercaptoethanol. 
Transferred  NOESY  (trNOESY)  experiments  (300  ms  mixing  time)  were  used  to  detect 
whether  binding  occurred  upon  addition  of  40  molar  excess  of  PhIP  (120  pM)  to  a  3  pM 
solution  of  ER  in  NMR  buffer.  TrNOESY  is  useful  in  detecting  binding  when  the  ligand  is 
in  fast  exchange  between  the  free  and  bound  states  (10"2-10~6  M  range  for  the 
dissociation  constant).  Adding  molar  excess  of  E2  to  the  solution  of  ER  and  PhIP  results 
in  precipitation.  This  is  most  likely  due  to  the  fact  that  E2  strongly  binds  ERa.  To 
determine  if  E2  could  displace  PhIP  from  the  binding  pocket,  0.25  molar  ratio  of  E2  (1 .5 
pM)  was  added  to  the  6  pM  ERa  protein-240  pM  PhIP  solution.  300  increments  were 
collected  in  ti,  each  with  128  or  256  scans  and  1024  complex  data  points  collected  in  t2. 
The  total  amount  of  DMSO  in  the  protein/ligand  samples  did  not  exceed  2%.  NMR  data 
was  processed  using  VNMR  software  (Varian  Inc.,  Palo  Alto  CA),  and  the  2D  frequency 
domain  matrices  were  analyzed  using  FELIX  (version  97,  Accelrys,  San  Diego,  CA). 


Cell  Culture 


MCF-7  cells  were  obtained  from  American  Type  Culture  Collection  (ATCC, 
Manassas,  VA)  and  were  grown  in  DMEM  with  5%  FBS,  1%  non-essential  amino  acids, 
10  pg/ml  insulin,  2  mM  L-glutamine,  and  1%  penicillin/streptomycin.  All  tissue  culture 
supplies,  with  the  exception  of  the  charcoal  stripped  fetal  bovine  serum,  were  obtained 
from  Invitrogen  (Carlsbad,  CA).  Cells  were  maintained  at  37°C  with  5%  humidity.  For 
proliferation  assays  and  estrogen  responsive  reporter  assays,  Opti-MEM  (phenol  red- 
free  medium)  was  used  with  5%  charcoal  stripped  fetal  bovine  serum  (Clontech,  Palo 
Alto  CA)  to  minimize  estrogen-like  activity  contributed  by  the  phenol  red  and  serum 
(Soto  et  al.,  1995;  Zacharewski,  1998). 

Cell  proliferation  assay 

For  cell  proliferation  assays  3  x  103  cells  were  seeded  into  each  well  of  a  96-well 
tissue  culture  plate  in  Opti-MEM  medium.  Twenty-four  hours  after  plating,  the  medium 
was  replaced  with  medium  supplemented  with  increasing  concentrations  of  E2  (Sigma, 
St.  Louis  MO),  PhIP  (Toronto  Research  Chemicals,  Downsview,  Ontario),  MelQx 
(Toronto  Research  Chemicals,  Downsview,  Ontario)  and  IFP,  treated  with  the 
appropriate  concentration  of  vehicle  or  left  untreated.  IFP  is  a  natural  product  isolated 
from  a  heated  mixture  of  creatine,  glutamic  acid  and  sucrose  (Pais  et  al.,  2000). 
Heterocyclic  amines  were  initially  dissolved  in  DMSO  and  then  serially  diluted  in 
medium.  DMSO  concentration  was  never  greater  than  0.5%.  Ethanol  was  used  to 


dissolve  E2.  The  effect  of  the  heterocyclic  amines  was  determined  by  comparing  ceil 
growth  to  the  appropriate  vehicle-treated  cells  and  untreated  cells. 

The  effect  of  test  compounds  on  cell  growth  was  evaluated  48  and  72  hours  after 
the  initiation  of  treatment.  Cell  growth  was  quantified  using  an  Aqueous  Non- 
Radioactive  Cell  Proliferation  Assay  (Promega,  Madison  Wl)  with  absorbance  of  the 
wells  measured  in  a  standard  multi-well  plate  reader  at  490nm.  Briefly,  this  colormetric 
assay  measures  the  bioreduction  of  a  tetrazolium  compound  (3-(4,5-dimethylthiazol-2- 
yl)-5-(3-carboxymethoxyphenyl)-2-(4-sulfophenyl)-2H-tetrazolium,  inner  salt;  MTS)  to  a 
formazan  product.  The  conversion  of  MTS  is  directly  proportional  to  the  number  of  living 
cells  in  culture.  Each  experiment  was  performed  at  least  three  individual  times  with  4  - 
6  replicates  per  experiment.  Effect  of  chemical  was  determined  by  comparing  cell 
growth  to  the  appropriate  vehicle-treated  cells  and  untreated  cells.  Error  is  presented  as 
the  standard  error  of  the  mean.  Statistical  significance  was  assessed  using  Wilcoxon 
Rank  Sum  Test. 


o  Estrogen  Responsive  Reporter  Assay 
Heterocyclic  amines  were  tested  for  estrogenic  and  anti-estrogenic  activity  using 
a  standard  estrogen  responsive  reporter  plasmid  containing  three  vitellogenin  estrogen 
responsive  elements  (EREs)  upstream  of  the  luciferase  reporter  gene  (EREx3-Luc). 
The  EREx3-Luc  was  a  kind  gift  of  Diane  Klotz,  NIEHS.  EREx3-Luc  was  co-transfected 
with  a  commercially  available  control  renilla  reporter  plasmid  (pRL-TK,  Promega, 
Madison  Wl)  to  control  for  transfection  efficiency  and  for  the  normalization  of  the  results 
The  reporter  plasmids  were  transfected  into  MCF-7  cells  (Klotz  et  al.,  1996;  Nodland  et 


al.,  1997;  Reel  et  al.,  1996;  Safe  et  al.,  1998},  which  constitutively  express  ERaand 
have  well-characterized  responses  to  estrogenic  stimuli.  The  ability  of  the  compounds  to 
induce  luciferase  reporter  gene  activity  was  compared  to  10'9  M  E2  and  appropriate 
negative  controls,  as  previously  described  (Charles  et  al.,  2002). 

MCF-7  cells  were  plated  in  96-well  plates,  designed  specifically  for  use  in  a 
luminometer  (Perkin  Elmer,  Wellesley,  MA)  at  a  density  of  1X105  cells/well  in  Opti-MEM 
medium.  The  cells  were  left  untransfected  or  were  co-transfected  with  200  ng  each  of 
luciferase  reporter  and  renilla  control  plasmids  using  lipofectAMiNE™  2000  (Invitrogen, 
Carlsbad,  CA)  optimized  for  96-well  plate  experiments  according  to  directions.  After  a  5- 
hour  transfection,  the  untransfected  cells  were  incubated  with  just  culture  medium  to 
quantify  background  levels.  The  transfected  cells  were  incubated  with  the  same  culture 
media  alone  or  supplemented  with  heterocyclic  amines,  E2  or  both.  Control  cells  were 
incubated  with  a  corresponding  amount  of  compound  solvent  at  the  highest 
concentration. 

The  cells  were  incubated  at  37 °C  for  18  hours  and  then  assayed  for  luciferase 
and  renilla  activity  using  a  dual  luciferase/renilla  reporter  assay  kit  (Promega,  Madison, 
Wl)  according  to  directions.  Each  experimental  condition  was  duplicated  six  to  eight 
times  in  the  96-well  plate  and  each  plate  was  read  in  the  luminometer  twice.  Luciferase 
activity  was  normalized  to  the  renilla  controls  and  compared  to  10"9  M  E2,  standard  in 
estrogen  responsive  luciferase  assays.  Estrogen  receptor-mediated  activation  was 
determined  by  comparisons  to  negative  and  E2  controls.  All  experiments  were  repeated 
more  than  3  times  and  all  data  were  normalized  to  E2  (10'9M)  activation.  For 
competition  assays,  HAs  (alone  or  in  combinations)  and/or  E2  were  added 


simultaneously.  Error  is  presented  as  the  standard  error  of  the  mean.  Statistical 
significance  was  assessed  using  Wilcoxon  Rank  Sum  Test. 

Results 

Docking  Studies 

Initial  docking  of  E2  to  the  ERa  LBD  as  a  control  calculation  reveal  the  ligand 
pose  as  seen  in  the  1GWR  crystal  structure,  with  the  initial  position  of  the  ligand  at  least 
15  A  outside  of  the  native  binding  cavity.  PhIP  and  A/^-hydroxy-PhlP  were  manually 
docked  in  the  same  orientation  as  E2  to  match  the  crystal  structure  (Brzozowski  et  al., 
1997),  where  the  exocyclic  amine  and  imidazole  ring  were  overlaid  on  the  E2  hydroxyl 
group  and  ‘A-ring’,  respectively.  These  initial  structures  started  the  docking  trials  of 
PhIP,  A/^hydroxy-PhlP,  IFP,  and  MelQx,  using  the  same  grid  of  the  protein  as  for  E2. 
Eleven  groupings  contain  99%  of  all  docking  poses  for  each  HA.  For  E2  and  PhIP,  the 
E2  binding  cavity  accounts  for  -98%  and  -89%  of  the  total  binding  probability, 
respectively.  PhIP  binds  favorably  into  the  E2  binding  site  although  the  calculated 
probability  of  binding  is  lower  when  compared  to  E2  (Fig.  2A).  Several  poses  of  PhIP 
are  observed  in  the  E2  binding  site.  One  pose  is  rotated  by  180°  with  the  phenyl  ring 
pointing  toward  residues  Glu353  and  Arg394.  However,  the  most  favorable  pose  places 
the  exocyclic  amine  proximal  to  the  hydrogen  bonds  between  Glu353  and  Arg394, 
similar  to  the  initial  docking  position  and  E2  pose  in  the  co-crystal  (Brzozowski  et  al., 
1997)  (Fig.  2B).  Side-chains  from  Leu387,  Leu391,  and  Phe404  make  hydrophobic 
contacts  with  each  face  of  the  imidazo  ring.  A  possible  hydrogen  bond  between  the 
exocyclic  primary  amine  and  the  carbonyl  of  Leu387  is  also  observed.  The  donor-H- 


acceptor  angle  is  greater  than  150°,  and  the  hydrogen  bond  distance  is  less  than  2.6  A 
(Fig.  2B).  The  phenyl  ring  of  PhIP  is  further  stabilized  by  hydrophobic  contacts  with 
Leu346  on  one  face  and  by  Met343  and  Met421  around  the  edge  of  the  ring.  The  face 
of  the  phenyl  ring  opposite  Leu346  is  open  to  solvent.  Furthermore,  the  phenyl  ring  is 
turned  relative  to  the  rest  of  the  molecule,  which  allows  for  a  better  fit  to  the  binding 
cavity  wall  and  is  consistent  with  the  geometry-optimized  structure  (Fig.  2B). 

Two  significant  groupings  of  //-hydroxy-PhIP,  IFP,  and  MelQx  are  observed,  one 
in  the  native  E2  binding  site  and  another  1 1  A  away,  close  to  helix  5.  The  lowest  energy 
pose  of  /'/-hydroxy- PhIP  is  in  the  E2  binding  site  and  is  in  the  same  orientation  as  PhIP 
(i.e.  all  hydrophobic  contacts  were  the  same).  However,  the  hydroxyl  group  does  not 
form  any  hydrogen  bonds  with  protein  side-chains  or  backbone.  Instead  of  forming  any 
obvious  hydrogen  bonds,  the  hydroxyl  group  was  further  buried  into  a  solvent  channel 
under  the  Glu353***Arg394  hydrogen  bond  bridge. 

The  planar  and  rigid  IFP  molecule  assumes  the  same  basic  orientation  of  PhIP 
and  //-hydroxy-PhIP;  the  exocyclic  amine  of  IFP  is  pointed  toward  the  Glu353»**Arg394 
hydrogen  bond  bridge.  A  hydrogen  bond  is  also  observed  between  the  primary  amine  of 
IFP  and  the  carbonyl  oxygen  of  Leu387.  The  binding  orientation  for  MelQx  in  the  E2 
binding  cavity  is  the  opposite  that  of  PhIP  and  //-hydroxy-PhIP  because  the  exocyclic 
amine  points  toward  His524,  similar  to  the  orientation  of  the  d-ring  of  E2  in  the  crystal 
structure  (1GWR).  We  observe  a  hydrogen  bond  with  the  side-chain  hydroxyl  group  of 
Thr347  and  the  exocyclic  amine  from  MelQx.  In  both  IFP  and  MelQx,  the  fused  rings  do 
not  appear  to  make  optimal  hydrophobic  interactions  in  comparison  with  PhIP  and  N2- 
hydroxy-PhlP. 


NMR  Ligand  Binding  Studies 

We  performed  transferred  NOESY  (trNOESY)  experiments  (Clore  et  a!.,  1982; 
Gronenborn  and  Clore,  1982;  Roberts,  1999)  to  detect  PhIP  binding  to  ERa  protein.  An 
observed  sign  flip  in  NOE  crosspeaks  of  PhIP  from  positive  (red)  to  negative  (black)  in 
the  TrNOESY  spectra  (Fig.  3A,  B),  is  an  indication  of  binding  to  ERa  protein.  Adding  E2 
to  the  PhIP/  ERa  mixture  causes  the  sign  flip  in  the  NOE  crosspeaks  to  revert  back  to 
positive  (red),  indicating  that  PhIP  is  no  longer  bound  to  ERa  (Fig.  3C). 

y 

MCF-7  Human  breast  cancer  cell  assays 

Incubating  ATCC  MCF-7  cells  with1x10‘9M  to1x10"7M  E2  for  48  hours  produces  a 
slight  increase  in  cell  proliferation  (7-12%)  that  is  significantly  different  from  vehicle- 
treated  cells.  By  72  hours  the  increase  in  cell  number  is  20%  above  control  cells  and  by 
96  hours  the  increase  has  risen  to  40%  (data  not  shown).  This  increase  in  cell  number 
is  consistent  with  other  reports  of  proliferation  response  for  this  MCF-7  cell  line 
(Rasmussen  and  Nielsen,  2002).  Adding  increasing  doses  of  PhIP  for  48  hours  also 
causes  a  7-10%  increase  in  proliferation  that  is  significantly  different  from  control  cells, 
but  not  significantly  different  from  E2  at  similar  concentrations  (Fig.  4).  By  72  hours  the 
increase  in  proliferation  reaches  20%  for  4x10'8M  PhIP,  which  is  not  significantly 
different  from  the  increase  in  cell  number  caused  by  a  similar  concentration  of  E2  (data 
not  shown)  cells  treated  with  A/^-hydroxy-PhlP  show  a  slight,  non-significant,  increase  in 
cell  growth  at  low  doses  (4x10'9M).  However  at  higher  concentrations  (4x10"6M),  cell 


growth  is  inhibited.  Treating  cells  with  concentrations  of  MelQx  and  IFP  higher  than 
4x10~8M  for  48  hours  causes  a  significant  decrease  in  cell  growth  (Fig.  4). 

A  standard  luciferase  reporter  assay  was  used  to  measure  the  effect  of  HAs  on 
ERa-dependent  transcription  activity.  This  assay  relies  on  the  ERa  normally  expressed 
by  the  MCF-7  cells.  Active  ER|3  is  not  expressed  in  MCF-7  cells  (Bardin  et  al.,  2004; 
Cestac  et  al.,  2005;  Paruthiyil  et  al.,  2004;  St-Laurent  et  al.,  2005),  and  therefore  ER(3 
interactions  cannot  contribute  to  the  results.  As  shown  in  figure  5A,  increasing  amounts 
of  PhIP  causes  a  dose-dependent  increase  in  ERa  activity  that  is  significantly  different 
than  untreated  cells.  In  contrast,  increasing  amounts  of  /'/-hydroxy- PhIP  and  the  other 
HAs  do  not  show  this  activation.  Competition  assays  were  used  to  determine  the  effect 
of  simultaneous  incubation  of  HAs  (4x10~7M)  and  10'9  M  E2  on  ERa  activation  (Fig.  5B). 
Treating  MCF-7  cells  with  PhIP  and  E2  produces  a  response  that  is  13%  less  than  E2 
alone.  //-Hydroxy-PhIP  inhibits  the  E2-mediated  activation  of  ERa  by  40%,  and  MelQx 
and  IFP  inhibit  ERa  activation  to  levels  that  are  significantly  lower  than  untreated  cells. 
Cell  viability  was  measured  under  identical  conditions  and  is  not  affected  by  incubation 
with  HAs  separately  or  in  combination  with  E2  at  these  concentrations  (data  not  shown). 

Effect  of  HA  combinations  on  ERa  activity 

Competition  assays  were  performed  with  PhIP,  //-hydroxy-PhIP,  and  MelQx  to 
determine  the  effect  of  simultaneous  exposure  to  4x10'7M  PhIP  and  increasing  amounts 
of  the  other  compounds  (Fig.  5C).  Treating  MCF-7  cells  with  either  a  1 :1  or  1 :3  molar 
ratio  of  PhIP  to  other  compounds  significantly  inhibits  ERa  transcription  activity  to  levels 
that  are  lower  than  cells  exposed  to  PhIP  alone. 


Discussion 


Heterocyclic  amines  are  known  to  be  potent  mutagens  and  rodent  carcinogens, 
but  the  complete  mechanism  of  carcinogenicity  for  these  compounds  has  not  been 
elucidated.  This  report  is  the  first  to  show  that  PhIP  binds  directly  to  the  ERa-ligand 
binding  domain  and  competes  with  estradiol  in  the  native  binding  cavity.  We  also 
demonstrate  the  novel  finding  that  A/^-hydroxy-PhlP,  the  PhIP  phase  I  metabolite,  and 
two  heterocyclic  amines,  MelQx  and  IFP,  inhibit  ERa  activation.  Among  the  HAs  tested, 
activation  of  the  receptor  is  limited  to  PhIP,  and  other  HAs  can  act  as  potential  anti¬ 
estrogens.  We  also  show  that  PhIP  increases  MCF-7  cell  proliferation  and  activates 
ERa  dependent  transcription,  in  agreement  with  results  published  by  Lauber  et  al. 
(Lauber  et  al.,  2004). 

PhIP  Binding  to  Era 

The  docking  of  E2  to  the  ERa  LBD  validates  our  model  of  ERa  LBD  and  the 
docking  algorithm  for  our  system  because  the  most  stable  E2  pose  found  was  the  same 
orientation  as  the  co-crystal  E2  conformation.  The  most  favorable  binding  pose  of  PhIP 
was  also  found  in  the  E2  binding  site,  even  though  the  probability  of  PhIP  binding  is 
lower  than  E2  binding  to  the  E2  site.  This  lower  binding  probability  is  consistent  with  the 
known  high  binding  affinity  of  E2  to  ERa  (Kd  ~  10'9-10'10  M)  and  our  NMR  results  show 
E2  competing  out  PhIP  (see  below). 


We  designed  an  NMR  competition  experiment  implementing  trNOESY 
techniques  that  provides  information  about  the  PhIP  binding  site  on  the  ERa  protein.  If 
PhIP  binds  the  same  site  as  E2,  then  adding  E2  to  the  PhIP/ERa  mixture  should 
displace  PhIP  since  E2  binds  more  strongly  to  ERa  than  PhIP  (as  predicted 
computationally).  If  PhIP  is  binding  in  a  different  site  than  E2,  then  the  negative  (black) 
PhIP-bound  peaks  should  remain.  When  E2  is  added  to  the  PhIP/  ERa  mixture,  a 
reversion  in  the  NOE  crosspeaks  from  the  negative  (black)  back  to  positive  (red)  is 
observed,  indicating  that  PhIP  is  no  longer  bound  to  ERa.  These  data  highly  suggest 
that  PhIP  is  binding  to  the  ligand-binding  domain.  However,  another  possibility  is  that 
PhIP  is  binding  in  another  site,  and  E2  is  preventing  PhIP  from  binding  because  of  a 
conformation  change  in  the  protein.  Although  this  possibility  may  exist,  the 
computational  prediction  of  PhIP  binding  in  the  E2  binding  site  and  our  in  vitro  data  on 
stimulated  cell  proliferation  (see  below)  leads  us  to  the  conclusion  that  PhIP  is  binding 
in  the  E2  binding  site  competitively. 

The  effect  of  the  HAs  on  the  proliferative  activity  of  whole  cells  was  tested  using 
a  modified  E-SCREEN  assay  (Rasmussen  and  Nielsen,  2002).  This  assay  is  based  on 
the  ability  of  estrogen-responsive  MCF-7  cells  to  proliferate  in  the  presence  of 
compounds  that  mimic  estrogen.  Although  a  standard  assay  for  screening  for 
estrogenicity,  the  proliferative  response  of  the  MCF-7  cells  varies  depending  on  the 
MCF-7  subline  employed.  In  these  studies  we  used  ATCC  wild-type  MCF-7  cells,  which 
have  been  shown  to  vary  in  their  response  to  estrogen  incubation  (Rasmussen  and 
Nielsen,  2002).  Using  these  cells,  we  found  that  micromolar  concentrations  of  PhIP 
were  able  to  stimulate  proliferation  up  to  40%  above  vehicle-treated  cells  over  a  72-hour 


incubation.  Lauber  et.  al.  determined  that  the  same  increase  in  cell  proliferation  could 
be  caused  by  only  nanomolar  concentrations  of  PhIP  (Lauber  et  al.,  2004).  These 
authors  used  MCF-7  cells  obtained  from  the  European  Collection  of  Cell  Cultures 
(ECACC)  in  their  assays;  presumably  this  MCF-7  cell  variant  is  more  sensitive  to 
estrogen  or  estrogen-like  compounds.  Nevertheless,  both  studies  demonstrated 
increased  cell  proliferation  when  incubated  with  PhIP. 

After  much  discussion  regarding  the  potential  estrogenicity  of  PhIP  and  its 
relevance  to  tumor  progression  (Gooderham  et  al.,  2002;  Lauber  et  al.,  2000;  Pfau  et 
al.,  2000),  we  and  others  (Lauber  et  al.,  2004)  now  confirm  that  PhIP  activates  ER  and 
stimulates  cell  growth.  Both  E2  and  PhIP  are  able  to  activate  ERa-dependent 
transcription  in  MCF-7  cells  transfected  with  a  standard  estrogen  responsive  reporter 
plasmid  containing  three  vitellogenin  EREs  upstream  of  the  luciferase  reporter  gene.  In 
other  words,  PhIP  is  able  to  use  the  endogenous  ERa  present  in  MCF-7  cells  to  effect 
transcription.  Thirty  minute  pre-incubation  with  the  complete  anti-estrogen  ICI  182,780 
inhibited  this  increased  activity  (data  not  shown)  (Lauber  et  al.,  2004),  further 
suggesting  that  PhIP  is  effecting  transcription  through  direct  interaction  with  ERa.  These 
results  and  the  fold-activation  of  ERa  transcription  activity  (compared  to  E2)  in  the 
presence  of  micromolar  concentrations  of  PhIP  are  in  agreement  with  previous  work 
(Lauber  et  al.,  2004). 

The  in  vitro  data,  without  the  supporting  NMR  and  computational  data,  could  also 
be  explained  by  PhIP  acting  on  the  ERa  pathway  at  points  upstream  or  downstream 
from  the  receptor  protein.  For  the  present  studies,  the  trNOESY  results,  together  with 
the  in  vitro  data  and  docking  results  strongly  support  the  idea  that  PhIP  and  E2  compete 


for  the  same  binding  cavity  within  the  ERa  LBD.  Taken  together,  these  results 
demonstrate  that  PhIP  initiates  an  estrogenic  response  in  target  cells  by  binding  directly 
with  the  ERa. 

A f-Hvdroxv-PhIP,  MelQx  and  IFP  Inhibition  of  ERa 

Close  examination  of  our  docking  data  shows  that  PhIP  makes  a  hydrogen  bond 
between  the  exocyclic-amine  and  the  carbonyl-oxygen  of  Leu387  (Fig.  2B).  A/^-Hydroxy- 
PhlP  shares  all  the  same  hydrophobic  contacts  as  PhIP.  However,  because  of  the  non- 
optimal  placement  of  the  hydroxy  group,  no  hydrogen  bonds  are  made  with  the  protein, 
decreasing  its  probability  of  binding  in  the  E2  cavity  (Fig.  2A).  These  data  suggest  that 
the  driving  force  in  binding  in  the  E2  cavity  is  hydrophobic  and  that  electrostatic 
interactions  such  as  hydrogen  bonds  play  only  a  secondary  role. 

In  the  cellular  assays,  N2-hydroxy-PhlP  inhibited  the  ability  of  E2  to  activate  ERa. 
When  compared  to  PhIP,  the  docking  data  predicts  that  the  binding  of  N2-hydroxy-PhlP 
should  be  weaker.  The  cell  assays  confirm  these  predictions.  When  N2-hydroxy-PhlP  is 
co-incubated  with  physiological  concentrations  of  E2,  much  of  the  ER  activation 
remains,  showing  that  N2-hydroxy-PhlP  is  not  a  particulary  effective  inhibitor. 

Like  PhIP  and  N2-hydroxy-PhlP,  MelQx  and  IFP  both  contain  an  exocyclic  amine 
with  a  methyl  group  at  the  adjacent  carbon.  However,  MelQx  and  IFP  contain  three 
fused-heterocyclic  rings,  which  limits  hydrogen  bonding  and  hydrophobic  interactions 
compared  to  PhIP.  The  low  barrier  rotation  of  the  phenyl  ring  in  PhIP  optimizes  non- 
bonded  contacts,  which  is  not  possible  in  MelQx  and  IFP.  Our  docking  data  suggest 
that  MelQx  and  IFP  do  not  bind  favorably  in  the  E2  binding  cavity  because  the  fixed 


conformation  of  the  fused-rings  does  not  allow  optimization  of  hydrophobic  contacts. 

The  lack  of  a  favorable  binding  conformation  is  confirmed  by  the  inability  of  these  to 
compounds  to  stimulate  reporter  gene  transcription  in  our  in  vitro  cell  assay  (Fig.  5). 

When  MelQx  and  IFP  are  co-incubated  with  physiological  concentrations  of  E2, 
ERa  activity  is  significantly  inhibited.  This  inhibition  is  identical  to  results  obtained  when 
the  compounds  are  added  alone.  Compared  to  N2-hydroxy-PhlP,  MelQx  and  IFP  are 
much  more  effective  inhibitors  so  that  E2  activity  remains  below  control  levels.  These 
results  suggest  that  the  mechanism  of  binding  for  MelQX  and  IFP  are  different  than 
PhIP  and  N2-hydroxy-PhlP;  MelQx  and  IFP  may  even  bind  in  other  regions  of  the  ERa 
protein. 

Implications  of  Diet  and  Metabolism 

The  estrogenic  potency  of  PhIP  is  comparable  to  other  environmental  estrogens 
found  in  the  diet.  Genistein,  a  known  phytoestrogen  found  in  soy  products,  was 
compared  to  PhIP  in  ERa  activation  assays  (data  not  shown).  At  low  concentrations 
(4X10'7M)  there  was  no  statistical  difference  in  the  activation  of  the  ERa  by  either 
compound,  but  higher  concentrations  of  genistein  elicited  a  much  greater  ERa  response 
than  higher  concentrations  of  PhIP.  When  investigating  compounds  from  food,  it  is 
difficult  to  predict  exposure  levels  for  individual  target  cells.  The  concentrations  of  PhIP 
examined  here  are  orders  of  magnitude  higher  than  what  a  single  cell  may  be  exposed 
to  after  a  meal  of  cooked  meat.  However,  chronic  exposure  to  PhIP  over  a  lifetime  may 
add  to  the  total  estrogenic  burden  of  the  body. 

In  addition,  it  is  not  clear  what  the  biological  consequence  of  being  exposed  to 
xeno-estrogens  may  be.  Although  exposure  to  environmental  estrogens  may  increase 


the  activity  of  endogenous  ERa,  exposure  to  other  compounds  may  inhibit  the  activation 
(anti-xeno-estrogens),  effectively  canceling  out  the  effect  of  the  xeno-estrogens.  In  fact, 
although  exposure  to  PhIP  may  be  stimulating  ERa  activity,  simultaneous  consumption 
of  the  other  HAs,  MelQx  and  IFP,  may  diminish  that  effect.  Our  previous  studies  have 
shown  that  when  meat  is  cooked,  the  ratios  of  formation  of  PhlP:MelQx  and  PhlP:IFP 
varies  according  to  cooking  conditions  (Knize  et  al.,  1998).  In  figure  5C,  1 :1  ratios  of 
PhlP:MelQx  and  PhlP:IFP  mixtures  significantly  inhibit  the  ability  of  PhIP  to  activate 
ERa.  If  the  ratios  are  increased  to  1 :3,  both  MelQx  and  IFP  completely  prevent  any 
stimulatory  action  by  PhIP.  When  cells  are  exposed  to  combinations  of  HAs  and 
metabolites,  as  would  happen  during  dietary  consumption  of  cooked  meat,  the  PhlP- 
mediated  activation  of  ERa  can  be  abolished.  Thus,  the  biological  consequences  of 
exposure  to  HAs  may  depend  upon  many  complex  factors,  including  the  ratio  of  HAs 
formed  in  the  meat,  how  well  each  compound  is  absorbed  from  the  digestive  tract,  and 
the  concentration  of  the  compounds  at  the  target  cells. 

Metabolism  of  the  compounds  will  also  contribute  to  the  complexity  of  the  overall 
human  exposure.  MCF-  7  cells  contain  active  P450  metabolism;  however  most  of  the 
activity  is  CYP1  Al .  Our  previous  investigations  have  shown  that  these  cells  do  not 
significantly  metabolize  PhIP  (K.  S.  Kulp,  unpublished  results),  suggesting  that 
metabolites  of  the  HAs  are  not  confounding  our  measurements  of  the  effects  of  the 
compounds  on  cell  proliferation  and  ERa  activity.  However,  in  humans  consuming 
cooked  meats,  hepatic  and  extra-hepatic  metabolism  of  PhIP  to  A^-hydroxy-PhlP  is 
extensive.  In  fact,  only  a  small  percentage  of  an  ingested  and  absorbed  dose  of  PhIP  is 
excreted  as  the  parent  compound;  the  rest  is  metabolites  (Malfatti  et  al.,  1999).  Our 


results  show  that  co-incubation  of  PhIP  and  AZ-hydroxy-PhlP  inhibits  the  stimulatory 
effect  of  PhIP  on  the  ERa  and  that  the  amount  of  inhibition  depends  upon  the  amount  of 
/V^-hydroxy-PhlP  present  in  the  incubation.  This  suggests  that  the  estrogenic  potential 
of  PhIP  may  depend  in  large  part  on  the  extent  of  Phase  I  activation  of  the  compound, 
either  in  the  liver  or  at  the  level  of  the  target  cell. 

Conclusions 

Here  we  show  by  experiment  and  computational  methods  that  the  food  mutagen 
PhIP  can  activate  human  ERa  by  binding  to  the  LBD  (Figs.  2,  3),  and  that  PhIP 
stimulates  MCF-7  breast  cancer  cell  proliferation,  suggesting  a  mechanism  for  tissue- 
specific  carcinogenesis  of  this  carcinogen.  We  also  show  that  IFP,  MelQx,  and  the 
primary  hydroxylated  metabolite  of  PhIP  do  not  stimulate  breast  cancer  cell  proliferation 
(Fig.  4).  Moreover,  MelQx,  IFP  and  A/^-hydroxy-PhlP  inhibit  estrogen  receptor  activation 
(Fig.  5).  Together,  these  results  imply  that  these  dietary  constituents  may  play  an  active 
role  in  both  activating  and  inhibiting  hormone  sensitive  cancers.  More  importantly, 
understanding  the  mechanisms  of  action  for  the  inhibitory  HAs  may  lead  to  potent 
therapeutics  against  breast  cancer. 
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Figure  Legend 

Figure  1 .  Structures  of  heterocyclic  amines  and  E2  (A)  focused  on  in  this  study  and  1 D- 
’H  spectra  of  5.95  mM  PhIP  (B)  and  132  mM  E2  (C)  dissolved  in  deuterated  DMSO 
(1 00%)  showing  the  relative  position  of  the  NMR  peaks  for  each  compound. 

Figure  2.  Boltzmann  probability  distributions  for  each  ligand  in  the  E2  binding  cavity  (A). 
Cut  away  view  of  the  ligand-binding  cavity  of  ERa  (B)  with  lowest  energy  docking  poses 
of  E2  (purple)  and  PhIP  (colored-sticks).  The  surface  is  colored  according  to  residue 
type;  nonpolar  (white),  polar  uncharged  (green).  Key  hydrogen  bonding  residues  are 
shown  in  colored-stick  representations  (Glu  353,  Arg  394,  Leu  387,  His  524).  Figure 
was  made  with  VMD  (Humphrey  et  al.,  1996)  and  rendered  with  Raster3d  (Merritt  and 
Bacon,  1997). 

Figure  3.  The  expanded  region  of  a  900  ms  mixing  time  NOESY  spectrum  (A)  of  5.95 
mM  PhIP  in  DMSO  (<2%)  exhibits  weak  positive  NOE  crosspeaks  (boxed  in  red).  The 
expanded  regions  of  300  ms  mixing  time  NOESY  spectra  of  ERa  with  (B)  addition  of  40 
fold  molar  excess  PhIP  indicating  that  PhIP  binds,  as  evidenced  by  the  flip  in  sign  to 
negative  crosspeaks  (boxed  in  black),  and  (C)  after  addition  of  0.25  molar  ratio  of  E2  to 
a  1 :40  ERa  +  PhIP  sample  showing  that  PhIP  does  not  bind  (red  peaks)  when  E2  is 
present. 

Figure  4.  Effect  of  PhIP,  Af-Hydroxy-PhIP,  MelQx,  IFP  and  E2  on  MCF-7  cell 
proliferation.  Cells  were  treated  with  compounds  for  48  hours.  Error  is  presented  as 
standard  error  of  the  mean.  *  =  significantly  different  from  control,  p<0.01 .  For  PhIP 
treatments  only,  *  =  significantly  different  from  control,  p<0.01  and  not  significantly 
different  from  estradiol  treatment  at  a  corresponding  molar  concentration.  X  Estradiol 
♦  PhIP,  ■  A/^-hydroxy-PhlP,  ▲  MelQx  and  •  IFP. 

Figure  5.  The  effect  of  heterocyclic  amines  on  estrogen  receptor  activation.  Cells  were 
treated  with  compounds  as  described  in  Methods.  Error  is  presented  as  standard  error 
of  the  mean.  Compound  concentrations  are  molar  and  E2  concentration  is  10"9M  in  all 
experiments.  Increasing  concentrations  of  Has  (A).  A/^-Hydroxy-  PhIP,  MelQx  and  IFP 
are  significantly  different  from  PhIP-treated  cells  at  the  same  concentrations  (p<0.01). 
All  compound  treatments  are  significantly  different  from  control  and  estradiol  (p<0.01). 
Competition  assays  with  estradiol  (B).  Cells  were  treated  with  4x10'7M  compound  alone 
or  compound  with  10'9M  E2  (hatched  bars).  Effect  of  combinations  of  heterocyclic 
amines  on  estrogen  receptor  activation  (C).  Simultaneous  addition  of  PhIP  and  N2- 
hydroxy-  PhIP,  or  PhIP  and  MelQx  are  significantly  different  from  PhIP  alone  (*) 
(p<0.01)  All  compound  treatments  are  significantly  different  from  control  and  estradiol 

(p<0.01). 
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Appendix  B 

Statement  of  Work 

Computational  Characterization  and  Prediction  of  Estrogen  Receptor  Coactivator 

Binding  Site  Inhibitors 

Task  1 :  Computationally  predict  the  relative  binding  of  CBIs  in  the  coactivator  site  and 
the  CBI  binding  properties.  (Months  1-12) 

(1 .1 )  A  minimum  of  three  CBIs  and  a  minimum  of  three  estrogenic  compounds  will  be 
docked  into  the  coactivator  site  and  the  estradiol  site,  respectively,  using  AutoDock. 
(Months  1-2) 

(1.2)  The  resulting  combination  of  CBIs  and  estrogenic  compounds  will  be  rank  ordered 
by  effective  binding.  (Months  3-4) 

(1 .3)  From  the  resulting  docked  compounds,  molecular  dynamics  (MD)  simulations  will 
be  run  on  the  energetically  favorably  bound  CBIs  to  investigate  their  effects  on  protein 
binding  and  stability.  MD  simulations  will  be  carried  out  using  NAMD  and  each 
simulation  will  be  run  for  2  ns.  (Months  4-1 1 ) 

(1.4)  The  resulting  combination  of  CBIs  and  estrogenic  compounds  will  be  re-ranked 
based  on  the  MD  simulations.  (Month  12) 

Task  2:  Biologically  assay  the  effect  of  the  CBI  on  the  estrogen  receptor  activation 
(Months  1-12) 

(2.1 )  Order  and  obtain  the  CBIs  and  estrogenic  compounds.  (Month  1 ) 

(2.2)  Measure  estrogen  receptor  activation  by  transfecting  a  reporter  plasmid  that 
contains  3  vitellogenin  of  estrogen  response  elements  upstream  of  a  luciferase  reporter 
gene  in  MCF-7  human  breast  cancer  cells.  Receptor  activity  will  be  measured  in  the 
presence  estradiol  and  each  individual  CBI.  (Month  2-3) 

(2.3)  Compare  the  effect  of  the  CBIs  on  estradiol  to  their  effect  on  the  other  estrogenic 
compounds  analyzed  in  the  docking  program.  (Months  4-6) 

Task  3:  Spectroscopically  measure  the  binding  of  the  CBI  to  the  ER  LBD. 

(3.1)  Obtain  the  compounds  used  in  the  biological  assay.  (Month  1) 

Order  and  obtain  a  pure  sample  of  the  human  ERa  (Month  1) 

(3.2)  If  the  CBIs  demonstrate  the  ability  to  modulate  estrogen  receptor  activation  in  the 
cells  assays  described  in  Task  2,  they  will  be  assayed  for  binding  to  Era  by  themselves, 
using  a  transferred  nuclear  Overhauser  effect  spectroscopy  (trNOESY)  assay.  (Months 
4-6) 

(3.3)  CBIs  that  bind  ER  and  demonstrate  activity  in  Task  2  will  be  assayed  against 
estradiol  in  a  competition  NMR  assay  to  determine  if  they  bind  in  the  presence  of 
estradiol.  If  they  bind  in  the  presence  of  estradiol,  then  the  CBIs  are  binding  to  another 
site.  (Months  6-8) 

(3.4)  To  rank  order  the  binding  of  CBIs,  CBIs  will  be  added  individually  to  ERa,  using  an 
NMR  assay  (Months  8-10) 

Task  4:  The  results  will  be  put  into  a  draft  of  a  paper  to  be  submitted  later  (Month  12) 


